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Collapse of a skyrmion due to the discreteness of a crystal lattice in isotropic two-dimensional 
ferro- and antiferromagnets has been studied analytically and by numerical solution of equations 
of motion for up to 2000x2000 classical spins on a square lattice coupled via Heisenberg exchange 
interaction. Excellent agreement between analytical and numerical results has been achieved. The 
lifetime of the skyrmion scales with its initial size, Ao, as (Ao/a) 5 in ferromagnets and as (Ao/a) 2 ' 15 
in antiferromagnets, with a being the lattice parameter. This makes antiferromagnetic skyrmions 
significantly shorter lived than ferromagnetic skyrmions. 
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Skyrmions [TJ |2] are tope-logically stable configurations 
of a fixed-length three-component vector field n(r) in the 
coordinate space of two dimensions. Due to the con- 
straint n 2 = 1 the n-field has two independent compo- 
nents. This permits unique mappings of n = (n x , n yi n z ) 
onto r = (x,y), described by classes of homotopy [3]. 
Each homotopy class corresponds to a non-trivial field 
configuration characterized by a conserved topological 
charge. The emergence of a conserved charge from a 
continuous field theory prompted numerous studies of 
skyrmions in problems of high-energy and condensed 
matter physics |4J. They include cosmology |5J, Bose- 
Enstein condensates [6], quantum Hall effect [7J [9] and 
anomalous Hall effect [8j, liquid crystals |10| . 

The interest to skyrmions in ordered spin systems 
had received much attention soon after the discovery 
of high-temperature superconductivity in copper oxides 
dHMHIMSj, and further explored recently [M22]. It 
is related to the fact that superconductivity in copper 
oxides occurs in doped Cu02 layers that, when undoped, 
are square lattices of antiferromagnetically ordered spins. 
Initially there was some hope that interaction of elec- 
trons and holes with skyrmions could play some role in 
Cooper pairing but this was never successfully demon- 
strated. Some indirect evidence of skyrmions in the mag- 
netoresistence of lanthanum copper oxide has been re- 
cently reported [23] but direct observation of skyrmions 
in 2d antiferromagnetic lattices is still lacking. 

In a continuous field model like, e.g., the non-linear 
a-model, the ground-state energy of the skyrmion does 
not depend on its size, A. This follows from the invari- 
ance of the model with respect to the scale transforma- 
tion r — >• &r, where k is an arbitrary constant. If the 
skyrmion lives on a lattice, however, the scale invariance 
becomes broken due to the presence of a lattice param- 
eter a. Thus the energy of the skyrmion depends on its 
size. This, in general, must lead to the collapse or expan- 
sion of the skyrmion, making it unstable. The nature of 
the exchange interaction on a lattice makes the skyrmion 
energy decreasing with its size, that leads to skyrmion 



collapse. A number of authors looked for interactions 
that could stabilize skyrmions in 2d ferromagnets [24U26] . 
It was argued that anisotropic crystal field added to the 
isotropic exchange model may, in principle, dynamically 
stabilize the skyrmion. In reality, however, anisotropic 
interactions are of relativistic origin, while the lattice ef- 
fect that leads to the collapse of the skyrmion is of the 
exchange origin and thus much greater. Therefore, it is 
important first to understand what is the mechanism of 
skyrmions collapse in a generic exchange model. 

In this Letter we study the dynamics of skyrmions 
and the dependence of their collapse time t c on their 
initial size in a 2d square lattice of classical spins cou- 
pled via Heisenberg ferromagnetic (FM) or antiferromag- 
netic (AFM) exchange interaction. The accuracy of the 
continuous approximation increases with the size of the 
skyrmion, A. One should, therefore, expect that the lat- 
tice skyrmion becomes stable in the limit of A — >> oo. We 
find that t c of the AFM skyrmion scales as t c oc (Ao/a) 215 
with its initial size Ao- We compute the dynamics of the 
collapse using both the analytical field model for the Neel 
vector and a direct numerical calculation on lattices of 
up to 2000x2000 exchange-coupled spins. The two ap- 
proaches show excellent agreement with each other. For a 
2d ferromagnet we obtain (up to logarithmic corrections) 
the (Ao/a) 5 scaling of the lifetime. This makes skyrmions 
significantly shorter lived in a 2d AFM than in a 2d FM. 

We begin with an antiferromagnet described by the 
Hamiltonian for the Neel vector L: 



Ho = l -JS 2 



dxdy 
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Here L is normalized as L 2 = 1, (VL) 2 = (d x L) 2 + 
(c^L) 2 , JS 2 > is the exchange energy associated with 
the interaction of spins of length S, and c is the speed of 
AFM spin waves that equals 2y/2 Ja/h in a square lattice. 
The term with L 2 can be understood as a kinetic energy 
responsible for the inertia of antiferromagnets. 

The absolute minimum of the energy corresponds to 
the uniform AFM background, L = const. Non-uniform 
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configurations of L are characterized by the topological 
charge 



(2) 



that takes values Q = 0, ±1, ±2, — Within, e.g., the 
homotopy class Q = — 1 the minimum energy, static con- 
figuration is a skyrmion given by 
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where r 
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- y 2 . Its energy, E = AttJS 2 , is indepen- 
dent of A. 

Equation ([!]) can be derived from the Heisenberg 
exchange interaction between nearest-neighbor classical 
spins |s A | = Is^l = 1, 
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where A and 5 denote AFM sublattices and ' = 
—Sl~L/S(Ss B,A ) = —S^2ij JijS^ ,A are the effective fields 
acting on the spins. As spins in each sublattice rotate 
smoothly through space, one can expand the effective 
fields as 



Hf = -JS 



4sf +a 2 V 2 sf 
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and similar for Hf . This allows one to go over to 
the continuum description in which there are two spin 
fields s A and s B . Switching to the magnetization M = 
(s A + s B ) /2 and the Neel vector L = (s A - s B ) /2, sat- 
isfying M 2 + L 2 = 1 and M • L = 0, with the help of 
equations of motion hs A,B = [s A,B x H A ' B ] one obtains 



H — Hq 
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JS a / dxdy 



(6) 



which differs from Eq. ([T]) by the second term due to the 
discreteness of the lattice. If the size of the skyrmion 
A is large compared to a, this term can be treated as a 
perturbation. Using the "rigid" skyrmion profile of Eq. 
([3]), one obtains the energy due to this term 
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that violates the scale invariance of the skyrmion. Eq. 
([7]) can be interpreted as a potential energy responsible 
for the skyrmion collapse. During the collapse it is trans- 
formed into the kinetic energy defined by the integral of 
L 2 = 4r 2 (r 2 + A 2 ) A 2 . With account of energy conser- 
vation, Eq. ([6| gives 
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Skyrmion collapse in the AFM model: 
analytical results 




Figure 1: Collapse of antiferromagnetic skyrmions as de- 
scribed by the numerical solution of Eq. ([8}. 



where Ao is the initial size of the skyrmion and r max has 
been introduced because of the logarithmic divergence of 
the integral in the kinetic energy. The natural choice is 
r max = Ao+ct, which describes a front of AFM spin waves 
propagating away from the collapsing skyrmion. This 
is confirmed by direct numerical calculations, see Fig. [4] 
below. The logarithmic terms with time-dependent r max 
require numerical integration of Eq. The resulting 
collapse curves are shown in Fig. [T] 

We now turn to the direct numerical solution of 
the dynamics of the skyrmion provided by the micro- 
scopic Hamiltonian Q. The dynamics is determined 
by the coupled equations of motion for spins, fast = 
— [&i x 6H/6(Ssi)]. We chose initial state as a staggered 
skyrmion texture, s st , given by Eq. (|3| for the A sublat- 
tice and by the same formula but with a minus sign for 
the B sublattice. The size of the skyrmion numerically 
can be defined as A^ = (m — l)(2 m 7r) _1 J2i (1 — ^li) 7 " \ 
where m > 1 is an integer. If one replaces summation by 
integration over dxdy /a 2 and uses the skyrmion texture 
(|3| for sf, this formula becomes an identity, A m = A. The 
results presented below have been obtained with m = 4. 
Other options make little difference. 

As the dynamics of the skyrmion is entirely due to 
small terms arising from the lattice discreteness, the 
time dependence is slow and sufficient accuracy can be 
achieved even for a large time step of integration. In- 
creasing the step is limited by stability rather then by 
required accuracy. The challenge of the numerical solu- 
tion is the 1/r decay of the skyrmion profile that requires 
rather big lattice sizes even for moderate values of A/a. 
Free or periodic boundary conditions introduce spurious 
A-dependent energies that compete with the small energy 
due to the lattice discreteness, leading to the expansion 
of the skyrmion instead of collapse. To make boundary 
conditions more resembling an infinite lattice, we have 
included the missing outside neighbors of the bound- 
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■ Skyrmion collapse in an antiferromagnet 
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Figure 2: Skyrmion collapse in an antiferromagnet. Whereas 
the skyrmion size A is decreasing continuously, the topological 
charge Q decays only during a short final stage of the collapse. 




100 200 300 400 500 600 



tJ/h 

Figure 3: Skyrmion collapse in an antiferromagnet for differ- 
ent initial skyrmion sizes. 



Figure 4: The front propagating from the center of antiferro- 
magnetic skyrmion at the beginning of its collapse. 
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Skyrmion collapse in the ferromagnetic model 
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Figure 5: Skyrmion collapse in a ferromagnet. The collapse 
time scales as t c oc Aq. 



ary spins with the values approximated by the second- 
order extrapolation from the inside of the working re- 
gion. Still, the lattice size has to be large: 1000x1000 
for X /a up to 16 and 2000 x 2000 for A /a = 18 and 
20. The program was implemented in Wolfram Mathe- 
matica with a compiled vectorized fixed-step fourth-order 
Runge-Kutta routine. One AFM skyrmion-collapse event 
required about one hour computer time. 

The collapse of an AFM skyrmion with Ao/a = 15 is 
shown in Fig. [2] Whereas the skyrmion size A is decreas- 
ing continuously, the topological invariant Q changes only 
during a short final stage of the collapse when the contin- 
uous approximation fails. Fig. [3] shows skyrmion collapse 
curves for different values of Ao/a. For Ao/a = 18 the lat- 
tice size of one million spins is too small and computation 
with four millions of spins is needed. For Ao/a = 16 these 
both lattice sizes yield the same collapse curve. These re- 
sults compare very well with the semi-analytical results 



shown in Fig. [T] The collapse time can be fitted as 
t c ex Aq' 15 in this range of Ao- The considerable deviation 
from the square law can be traced back to the logarith- 
mic term in Eq. Fig. [4] shows \dL/dt\ in an anti- 
ferromagnetic skyrmion at different times. The region of 
skyrmion motion where \dL/dt\ > is expanding with 
the speed of antiferromagnetic spin waves c. The reason 
for this is that the lattice-discreteness terms that drive 
the skyrmion collapse are very short ranged while the 
skyrmion itself is long-ranged. The action of the former 
is transferred to the whole skyrmion with a speed c in 
accordance with the causality. The front position can be 
estimated as r max = Ao + ct, as was argued after Eq. (|8|. 

Along the same lines we have numerically studied the 
dynamics of ferromagnetic skyrmions. It turns out to be 
much slower than the collapse of AFM skyrmions, so that 
up to one day of computations is needed for one collapse 
event. Fig. [5] shows time dependences of the size of FM 
skyrmions during the collapse. The collapse time scales 
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The Aq scaling of the collapse time of the FM skyrmion 
can be qualitatively understood as follows. The exchange 
interaction conserves the total spin of the system. The 
infinitesimal increase of the (negative) skyrmion spin in 
the course of its collapse is 



dS = S 



J a* 



ds z 
~dX 



dX 



o „XdX . R 
-8tt5— 5- In T , 
a z X 



dX > 0. 



(9) 

Here we used s z in the skyrmion form given by Eq. ([3| 
and introduced the long range cut-off R. Because of 
the conservation of the total spin, the increase of the 
skyrmion spin by dS generates dS magnons. Since in 
this process the spin is being carried by large distances, 
the skyrmion collapse is very slow. The average energy 
of emitted magnons can be estimated as Huj ~ —hX/a. 
This yields the emitted magnon power 
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fko— = 8nhS ^-Ln- 



dt 



A 



(10) 



On the other hand, the rate of change of the energy ^ 
due to discreteness of the lattice is £discr oc A. From the 
energy conservation, £discr + P = 0, one obtains 
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yielding the collapse time 
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The condition Huj <C SJ for the energy of the magnons 
translates to 5(A/a) 4 \n(R/X) ^> 1, which is well satisfied 
during the collapse. 

In conclusion, we have studied the collapse of 
skyrmions due to the discreteness of the lattice in generic 
models of isotropic 2d ferro- and antiferromagnets with 
Heisenberg exchange interaction. The results obtained 
within continuous field model are in excellent agreement 
with the direct numerical calculation on lattices of up to 
2000x2000 coupled spins. The collapse time of antiferro- 
magnetic skyrmions obtained by both methods scales as 
(Ao/a) 215 . For ferromagnetic skyrmions, the numerical 
calculation gives the (Ao/a) 5 scaling of the collapse time. 
It is explained by the emission of magnons. Thus, AFM 
skyrmions are much shorter lived than FM skyrmions. 
This can be understood in the following terms. The col- 
lapse of an AFM skyrmion occurs via transformation of 
its potential energy due to the discreteness of the lattice 
into the kinetic energy defined by L 2 . The FM skyrmion 
does not possess such a kinetic energy, so that its po- 
tential energy has to be dissipated into magnons, which 
is a much slower process. In the expression for t c the 



time constant in front of the power of the ratio Xq/cl is 
of order H/(JS). For, e.g., JS ~ 100 K and A - 10a, 
this gives t c ~ 10ns for the lifetime of the skyrmion in a 
ferromagnet and t c ~ lOps in an antiferromagnet. 
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